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Abstract 

This is a survey of results about permanental processes, real valued 
positive processes which are a generalization of squares of Gaussian pro- 
cesses. In a certain sense the symmetric positive definite function that 
determines a Gaussian process is replaced by a function that is not nec- 
essarily symmetric nor positive definite, but that nevertheless determines 
a stochastic process. This is a new avenue of research with very many 
open problems. 

1 Introduction 

Let T be an index set and G = {G(x);x € T} be a mean zero Gaussian 
process with covariance U = {U(x, y);x,y € T}. The covariance U completely 
determines G. If instead of G we decided to consider its square, i.e., G 2 := 
{G 2 (x); x € T}, we would not find this difficult, since we could simply obtain 
information about G 2 from our extensive understanding of the properties of 
G. However, let us suppose that we were completely unaware of the theory of 
Gaussian processes. We then might begin our study of G 2 by considering its 
finite joint distributions, as given by their Laplace transforms, i.e., 

(-^E^ 2 ^))) = | / + I r[ i/ 2 

for all x\, . . . , x n in T, where / is the n x n identity matrix, a is the diagonal 
matrix with (a^j = a^), € R+ and T = {T(xi,xj)} is an n x n matrix, 
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that is symmetric and positive definite. (Of course, in truth, we do know 
the theory of Gaussian processes so we know that {r(xj, xj)} is actually the 
covariance of the mean zero Gaussian vector, (G(xi), . . . ,G(x n )). That is, 
T = {U(xiXj)}fj = -±.) Nevertheless let us continue to pretend that we do not 
know the theory of Gaussian processes. A natural problem for us to consider 
is to find necessary and sufficient conditions for G 2 to have continuous paths. 
Considering that the only unknown on the right-hand side of (jl.ip is the matrix 
r, the answer would have to be given in terms of V. 

Lets now suppose that we have obtained necessary and sufficient conditions 
for G 2 to have continuous paths. We might then ask whether we could relax 
the conditions that T is symmetric and positive definite in (jl.ip and that the 
power of the determinant is 1/2. For example we might ask when we can 
define a real valued positive stochastic process 9 = {0 x ,x € T} with finite 
joint distributions that satisfy 

£ (»e(-!g aA -))=_W' (1.2) 

where T = {T(xi, xj)}fj =1 is an n x n matrix and j3 > 0. Of course, aside 
from dimension 1 where 9 Xl has a gamma distribution, for f3 > 0, it is not 
at all clear whether (|1.2p makes sense, i.e., the right hand side may not be a 
Laplace transform. 

In 1997, D. Vere- Jones, [13] , initiated the study of random vectors with 
Laplace transform given by (11.2p . He also briefly considered stochastic pro- 
cesses for which (II. 2D gives finite joint distributions. We call these random 
vectors and stochastic processes /3-permanental processes. 

In [13] Vere-Jones gives necessary and sufficient conditions in terms of T 
and a for (|1.2p to exist, i.e., that the right-hand side is a Laplace transform. 
Unfortunately, at least to us, it seems impossible to verify them, except in a 
very important case that was observed by N. Eisenbaum and H. Kaspi [5j. This 
is the case when T is the potential density of a transient Markov process on T, 
or, equivalently, when all the finite dimensional matrices T are M-matrices, 
or, equivalently, when all finite dimensional vectors (9 X1 , . . . , 6 Xn ) are infinitely 
divisible. None of these three mathematical concepts is generally considered 
by experts in Gaussian processes. One purpose of this review is to convince 
the reader that they should be. 

Several researchers have been looking for something like permanental pro- 
cesses for a long time. The Dynkin Isomorphism Theorem, (DIT) relates the 
squares of Gaussian processes to the local times of Markov processes, which 
necessarily have symmetric 0-potential densities. An outstanding problem 
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was to find a version of Dynkin's theorem when the O-potential densities of 
the Markov processes are not symmetric. This is no doubt what really inspired 
Eisenbaum and Kaspi. They give a generalization of the DIT, or to be more 
precise, of a version of the DIT given in [6], in [5j Corollary 3.5]. Of course, 
the generalization of the process that is Gaussian squares, is the permanental 
process. The DIT enables one to use results about Gaussian processes to ob- 
tain results about the local times of associated Markov processes. This is what 
[TT] is devoted to. In [121 Theorem 1.3] Theorem 13.11 in this paper is used to 
obtain a sufficient condition for the continuity of local times of Markov pro- 
cesses without requiring that they have symmetric O-potential densities. This 
result with a slightly more restrictive hypothesis was obtained in [4] by a dif- 
ferent method. We do not discuss Dynkin type isomorphism theorems further 
in this paper; our focus is on permanental processes. However, we advise any 
researcher who decides to work on permanental processes to study them. 

Another fascinating property of some permanental processes is that they 
can be represented as the local times of loop soups associated with Markov 
processes. This was shown for permanental processes associated with sym- 
metric Markov processes in [8], (see also [9]), and extended to permanental 
processes associated with more general Markov processes in |12| . 

In Section [2] we give Verre- Jones's definition of f3 permanents and his nec- 
essary and sufficient condition for the existence of /3-permanental vectors. In 
Section El under the assumption that (]1.2p defines a permanental process, we 
give sufficient conditions for them to be bounded or continuous. In Section|3]we 
summarize some of the results in [5] that show that /3-permanental processes 
exist if their kernels are the zero potentials of transient Markov processes. (In 
which case the permanental processes are infinitely divisible.) Levy processes 
killed the first time they hit zero, or at the end of and independent exponen- 
tial time, have zero potentials that are kernels of permanental processes. In 
Section [5] we obtain these potentials of Levy processes and thus get an idea of 
what the kernels of permanental processes, that are not the squares of Gaus- 
sian processes, look like. To be specific, in Section 15.11 we give the kernels 
of permanental processes that are generalizations of the squares of fractional 
Brownian motion. The study of permanental processes is just beginning. In 
Section [6] we list a few problems that interest us. 

This paper is based on Eisenbaum and Kaspi's seminal paper [5] and on 
|12j . It is primarily a survey except that Remark 13.11 Lemmas I4.2H4.4I and 
Theorem 14.21 are new results. Also the process FBMQ ' 3 appears for the first 
time in this paper. 
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2 Existence of permanental processes 



The square of a Gaussian process is a permanental process so there are plenty 
of examples. But we are interested here primarily in those that are not squares 
of Gaussian processes. To begin it is important to note that the right-hand 
side of (jl.2p is not unique with respect to T. If D is any diagonal matrix with 
non-zero entries we have 

| J + a T\ = \I+ uD- 1 TD\ = \I + aD- l Y T D\. (2.1) 

For a very large class of irreducible matrices T, it is known that these are 
the only sources of non-uniqueness; see [10J. On the other hand, in certain 
extreme cases, for example, if T± and T2 are n x n matrices with the same 
diagonal elements and all zeros below the diagonal, then |J+ari| = IJ + cJ^I- 
Therefore, a process that is the square of a Gaussian process can have a non- 
symmetric kernel. Nevertheless, as we show, there are many permanental 
processes that are not squares of Gaussian processes. 

The first step in the definition of permanental processes is the definition of 
the /3-permanent of a n x n matrix B = {6ij}" J=1 , which we denote by \B\p. 
It is 

\B\p = Y J P m{<T) K^2,n---K,i n , (2.2) 
where the summation is taken over all distinct permutations 

.= (. 1 ' 2 "-") (2.3) 

of the indices 1, 2, . . . , n and m(a) denotes the number of cycles into which a 
can be decomposed. 

The name /3-permanent is taken from the fact that \B\i is the permanent 
of B. 

For /3 > we define what it means for an n x n matrix B to be j3 positive 
definite. This is a rather complicated definition. To begin consider all multi- 
indices k = (k\, . . . , k n ), (where the ki are integers), and set |k| = k± + - • - + k n . 
For each k, let B(k) denote the |k| x |k| matrix with entries 

{B(k)} - = B pup . (2.4) 

where pi = 1 if 1 < i < k\, and pi = I if k\^\ < I < ki, if I = 2, . . . , n. The 
matrix B to said to be f3 positive definite if for all multi-indices |k| 

|B(k)|^>0. (2.5) 
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Clearly, to verify whether B is f3 positive definite one must check (|2,5p for an 
infinite number of matrices. 

In [13\ Proposition 4.5] Verre- Jones gives necessary and sufficient condi- 
tions for (jl.2p to be the Laplace transform of the vector (6 Xl , . . . , 9 n ) in terms 
of the modified resolvent matrix 

T r :=T(I + rT)- 1 (2.6) 

where r > 0, and T is the matrix in (jl.2p . Note that when T _1 exists, we have 

r r = (r- 1 + ri)- 1 . 

The next proposition is Proposition 4.5 in |13j . 

Proposition 2.1 For il.2\) to represent the Laplace transform of a non-negative 
random vector it is necessary and sufficient that for all r > 

(i) T r exists and is /3 -positive definite. 

(ii) det(J + rr) > 0. 

(Item (ii) is equivalent to: All the real, non-zero, eigenvalues of V are 
positive.) 

Verifying Proposition 12. II (i), requires verifying (|2.4p for an infinite number 
of matrices, which seems impossible unless all the entries of the matrix T r are 
greater than or equal to zero. In [5, Theorem 3.1] Eisenbaum and Kaspi 
point out that this is is the case for an important class of kernels T and that 
Proposition 12.11 (ii), also holds for these kernels. We take this up in Section 
01 In the next section we assume that a permanental process exists and give 
sufficient conditions for it to be continuous. 



3 Continuity of permanental processes 

A key observation that allows us to obtain a 'best possible' sufficient condition 
for the continuity of permanental processes is that the bivariate marginals of 
a 1/2-permanental process are squares of bivariate normal random variables. 
We proceed to explain this. 

For n = 2, (jl.ip takes the form 

E(^(~( ai e x + a2e v )\\ (3.1) 

= | /+ * r ji/ 2 = ( 1 + a i T ( x , x ) + "2r(y, y) 

+aia 2 {r(x,x)T(y,y) - T(x,y)T(y,x))y 1/2 . 
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An obvious necessary condition for this to exist is that |/+aT| > 0. Therefore, 
if we set «2 = in (|3.ip and take a\ sufficiently large we see that for any x G T 



T(x,x) > 0, 
and taking a\ = «2 sufficiently large, that 

r(x,x)T(y,y) - F(x,y)T(y,x) > 0. 



(3.2) 



(3.3) 



In addition, by [131 p. 135, last line], a necessary condition for the existence 
of any /3-permanental process is that for any pair x, y G T 



F(x,y)T(y,x) > 0. 

This allows us to define 

d(x, y) = 4^2/3 (t(x, x) + F{y, y) - 2 (T{x, y)T(y, x)) 1 ' 2 
It follows from (|3.2p - (|3.4p that for any pair x, y G T, the matrix 



1/2 



r(x,x) (r(x,y)r(y,x)) 1/2 ' 
(r(x,y)r(y,x)) 1/2 F(y,y) 



(3.4) 



(3.5) 



(3.6) 



is positive definite. Therefore we can construct a mean zero Gaussian vector 
{G(x),G(y)} with covariance matrix 



Note that 



E(G(x)G(y)) = (T(x,y)T(y,x)) 1/2 . 

2 .l/2 v^/2. 



{E(G(x) - G(y)) 



d(x,y), 



(3.7) 



(3.8) 



By p.ip the Laplace transform of {9 x ,0y} is the same as the Laplace 
transform of {G 2 (x) , G 2 (y)} . This gives the following critical lemma: 

Lemma 3.1 Suppose that 9 := {0 x ,x G T} is a 1/2-permanental process with 
kernel T as given in M.l}) . Then for any pair x,y, 



{6 x ,e y }±{G 2 (x),G 2 (y)} 



(3.9) 



where {G(x),G(y)} is a mean zero Gaussian random variable with covariance 
matrix given by \3. 
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For p > 1, let tp p {x) = exp(x p ) - 1 and L^ P ($7,J 7 , P) denote the set of 
random variables ^ : — > i? 1 such that (l£l/ c ) < oo for some c > 0. 
L^ p (fi, J 7 , P) is a Banach space with norm given by 

||CIU P =inf{ C >0:^ P (!e|/c) <1}. (3.10) 

We shall only be concerned with the cases p = 1 and 2. 
It follows from Lemma |3. II that 

110s - 9 y \\^ = \\G 2 (x) - G 2 (y)\\^ < CsnpT 1 / 2 (x,x)d(x,y) (3.11) 

This is a good start but not the result we want since we know that the con- 
tinuity condition when the permanental process is the square of a Gaussian 
process is the sufficient condition for continuity of the Gaussian process itself. 
This necessitates working with the Orlicz space L^ 2 . 

Unfortunately, \\G 2 (x) — G 2 (y)\\^ 2 = oo. We get around this by using a 
truncation argument introduced by Martin Barlow to obtain a sufficient con- 
dition for the continuity of local times of Levy processes [H [2] . The following 
lemma is [121 Lemma 3.3]. 

Lemma 3.2 Let 9 := {9 x ,x € T} be a 1/2-permanental process with kernel 
T. Then for all x,y € T and < A < oo 



x A A By A A 



A 1 / 2 A 1 / 2 



<d(x,y). (3.12) 

W2 



The proof of Theorem 13.11 below uses (|3. 12|) and standard ideas used to 
prove the continuity of Gaussian processes along with some of the ideas in 

mm- 

Let (T, p) be a separable metric or pseudometric space. Let B p (t, u) denote 
the closed ball in (T, p) with radius u and center t. For any probability measure 
H on (T, p) we define 



teT Jo 



1 



JT, P A a ) = sup / log — r- du. (3.13) 



n{BJt,u)) 



1/2 



We occasionally omit some of the subscripts T, p or /i, if they are clear from 
the context. 

In general, d(x, y) is not a metric or pseudometric on T. Nevertheless, we 
can still define the sets Bd(s,u) = {t £ T\d(s,t) < u}. We can then define 
JT,d,fi(o-) as in (|3. 13|) . for any probability measure p on B(T,d), the cr-algebra 
generated by the sets Bd(s,u). 
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Theorem 3.1 Let 9 = {9 X : x € T} be a 1/2-permanental process with 
kernelT satisfying sup x£T T(x,x) < oo. Let D denote the d-diameter ofT and 
assume that T is separable for d, and that there exists a probability measure \i 
on B(T, d) such that 

J d {D) < oo. (3.14) 
Then there exists a version 6' = {9' x ,x € T} of 9 such that for any xq € T 



sup 0^11^ < 4\\9' X0 \\^+C (supT{x,x)) jj(D), 

x&T \xeT J 



(3.15) 



where C is a constant. 
If 



lim J d {5) = 0, (3.16) 
5^0 



lim Ml = oo, (3.18) 



there exists a version 9' = {9 x ,x € T} of 9 such that 

lim sup \e' s (cj) -6' t {u)\ = 0, a.s. (3.17) 

d(s,t)<8 

If holds and 

s'^b 5 
then 

\9' — 9'\ ( \ 1//2 

lim sup — i-f r^— <30 sup0' a.s. (3.19) 

5->o J d (d(s,t)/2) ~ Vxgt J 

d(s,t)<8 

Continuity is not mentioned in Theorem 13.11 because we do not know 
whether d is a metric on T. If d(x, y) is continuous on T x T and there 
exists a probability measure /ionT such that (|3.16p holds then there exists a 
version 9' = {9' x : x € T} of 9 that is continuous almost surely. In this case 
(|3.19p gives a bound for the uniform modulus of continuity. It s not unreason- 
able to assume that d(x, y) is continuous on T x T. We show in [121 Lemma 
3.2] that when 9 is continuous on T almost surely, d{x, y) is continuous on 
T x T. 

1 /2 

It is easy to see why the term sup xgT (9' x ) ' appears on the left-hand 
side of (13. 19H . Consider a Gaussian process {G x : x € T} with modulus of 
continuity u(d(s,t)). (Suppose that (T,d) is a metric space.) Then 

— \q — Q.\ ( X 1 / 2 
lim sup —. s — ^ < lim sup ^ \2supG 2 x ) a.s. (3.20) 



<5^o SjteT u(d(s,t)) s->o SjteT u}(d(s,t)) \ .,._/■ 

d(s,t)<8 d(s,t)<8 
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(With mild regularity conditions the inequality can be replaced by equality.) 

Under the hypotheses of Theorem 13.11 an essentially optimal local modulus 
is also given for 9 = {9 X : x G T} in |12|. Theorem 4.2]. 

We say that a metric or pseudometric d\ dominates d on T if 

d(x, y) < y), Vx, y G T. (3.21) 

In the \12\ Section 5] we give several natural metrics that dominate d. For 
example 

\ inf T 1 / 2 (x,x)d(x,y) < ( E (0 X - 9 y f V 2 < ^supT 1 / 2 (x,x)d(x,y). (3.22) 



i / v ' v v — •*" / i ,. , 

If a metric di dominates d and is such that (T, <ii ) is separable and has 
finite diameter D then the results of Theorem 13.11 hold with d replaced by d\ 
and the version 9' is uniformly continuous on (T, d\). (See |X2(, Corollary 1.2].) 

Remark 3.1 We expect that Theorem 13.11 also holds for all /3-permanental 
processes, j3 > 0, but we can't show this. What we can show is that all 
the statements in Theorem 13.11 and in the discussion following it hold for all 
/3-permanental processes, j3 > 0, if Jt,p,h{o) in f)3. 13 j) is replaced by 

X T , p » = sup / flog ^j^y ) du. (3.23) 



teT Jo 



This follows from the next lemma and corollary used in conjunction with the 
proof of |12l Theorem 1.1]. 

Lemma 3.3 Let 9 := {9 x ,x G T} be a f3-permanental process, f3 < 1/2 with 
kernel T. Then for all x,y G T and < A < oo 

< 16(1 ~ P) fsu P r(x,x)) d(x,y). (3.24) 
VI o \xeT J 

Proof Let T(x, y) be the kernel of the 2-dimensional vector 9 = (9 X , 9 y ). Let 
G = (G x , Gy) be the 2-dimensional Gaussian vector with covariance given by 
()3.7p . By (13. ip we can assume that the kernel of 9 is given by the matrix in 
(|3,6p which we denote by T. It is well known that 

E (exp \-~{a x e s + aye v )Y\ = J7^f> ( 3 - 25 ) 
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is the Laplace transform of vector in for all f3 > 0. (See Lemma 14.11 and 
Remark l4.1i ) Denote this vector by {0 x ^,6 y ^). For (3 < 1/2 we have 

(9x,/3 ~ 0y,p) + (#2,1/2-/3 ~ Q'y,l/2-p) = G l ~ G y ( 3 - 26 ) 

where 1/2-^^1/2-/3) * s independent of (0 x ,/3,Gy,p), since (G x ,G y ) is the 
vector obtained from (|3,25|) when (3 = 1/2. Let EqEqi denote expectation 
on the product space of {0 X) p,9 y ^) x [0' x 1 / 2 _ / g, 0' y 1/2-/3) an< ^ ^ II ' Hp = 
{E e E e ,\ ■ \p) 1/p , p> 1. We have 

II - Q y ,p) + (^,1/2-/3 ~ ^,1/2-/3) Hp ( 3 - 27 ) 

i/p 

i/p 

■/3 - y y,l/2-/3)rj 
1' Ml 



EeEe>\{9x,p - y ,p) + (0'x,l/2-p ~ 0> y,i/2-p)\ P 

> [Eg\(e Xj p - y> p) + Egf(e' x>1/2 _p - 6*^,1/2-^)1^ 

> ll^x,/3 - + E(0' Xt i/2-p ~ v y,l/2-p)\\V 

> Px,P ~ Qy,p\\p ~ \ E e'( e 'x,l/2-P ~ Qy,l/Z-p)V 

Consequently 

\\0x,p ~ Gy,p\\p - W G l ~ G yWp + 1^(^,1/2-/3 ~ ^,1/2-/3)!- (3.28) 
Using (|3.25p we see that for any /3 < 1/2 

1^(^,1/2-/3-^,1/2-/3)1 = (l-2P)\T(x,x)-T(y,y)\ (3.29) 

< (1-20)||G£-GJ|| p , 

since for integers p > 1, 

(£|G* - G*| p ) 1/p > (IE (G* - GJ) \ P ) 1/P = \T(x,x) - T(y,y)\. (3.30) 
Therefore 

E(0 Xtf) - 6 y ^r < (2(1 - P)fE{G 2 x - G 2 y y, (3.31) 
which implies that 

110* < 2(1-/9)||GS-G?||^ (3.32) 

which gives (I3.24p . □ 
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Corollary 3.1 Let 9 := {9 x ,x G T} &e a /3-permanental process, (3 > 1/2, 
with kernel V . Then for some integer k > 1 

<l^( supr( „)W (3.33) 

Proof Set P = k/3' for some (3' < 1/2. Then since 

\\6x,p — Qy,p\\ipi — k\\d X) p> - Oy^'W^, (3.34) 

(|3T33|) follows from ([3T2i]l . □ 

In general, given the existence of a 1/2-permanental process with kernel 
r, we do not know whether there exists any other /3-permanental process with 
kernel T. If one does exist for rational ft, the next theorem shows that it 
inherits continuity properties of the 1/2-permanental process. 

We use the following inequality in the proof: 

Lemma 3.4 [3, Theorem 1.1.5] Let X\, . . . , X n be independent identically dis- 
tributed random variables in a normed linear space. Then for all 1 < k < n 
and all t > 



(3.35) 



P(\\X 1 \\>t)<9P[\\^Tx j \\>^\. (3.36) 

3=1 

Theorem 3.2 Let 9 = {9 X : x E T} be a 1/2-permanental process with kernel 
r that has continuous paths on (T, d) almost surely. Let 9p = {6 X r : x E T} 
be a /3-permanental process with kernel T for any rational (3 > 0. Then 9r has 
continuous paths on (T, d) almost surely. 

Proof Let T' be a finite subset of T and |j • || := sup^^)^ | • |. Suppose 

x,yeT' 

that f3 = p/q, for integers p and q. Then 




In particular 

I n 

PIWX-, II ->f\ < QP I II V XA\ > 

30 



q 2p 

E " e v^ j) ) II = HEW* - W» W H- (3-37) 

3=1 3=1 
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where the notation (j) denotes independent copies. By (j3.36j) and ()3.37|) 

(3.38) 



P (\\0x,p - v ,p\\ >t) 

2p 



< 



w HE' 



3=1 



<18PP(||^ V2 -V/ 2 II>6^) 



18pP 



( 

SUP 1^,1/2-^,1/21 > ~T 

. d(x,y)<6 UU F 



Since 6^1/2 has continuous paths it follows from (|3.38p that for any finite 
subset T' of T and e > we can take 5 sufficiently small so that 



Pi 



> e < e. 



(3.39) 



This implies that Op has continuous paths on (T, d) almost surely. 



□ 



4 Examples of permanental vectors and processes 

The matrix A is said to be an M matrix if 

(1) a i,j < for all i 7^ j. 

(2) A is nonsingular and A^ 1 > 0. 

Lemma 4.1 [5, Lemma 4-2] Let {Tjj} be a real positive non-singular n x n 
matrix. There exists a positive infinitely divisible random vector (9i, . . . ,9 n ) 
such that for every (ati, . . . , a n ) in i2" 

("5 1 -«■)) =^r^ <4I) 

i/ and on/y i/T -1 is an M -matrix. 

Remark 4.1 One can replace the 1/2 in (14. If) by any j3 > 0. To see this 
denote the right hand side of (|4.ip by F(a). A necessary and sufficient condi- 
tion for a function F(a) to be the Laplace transform of an infinitely divisible 
random variable is that logF(5), where S is a diagonal matrix with entries 
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(1 — Sj)i, t > and < Sj < 1, i = 1, . . . , n, has a power series expansion in 
(si, . . . , s n ) with positive coefficients for all £ sufficiently large. (See, e.g., [ITJ 
Lemma 13.2.2].) In (|3~Tjl 



Clearly, if this has a power series expansion in (si, . . . , s n ) with positive coef- 
ficients for all t sufficiently large, it does if 1/2 is replaced by any j3 > 0. 

Let X = {Xt,t £ be a transient Markov process with state space E 
and potential density u(x,y). 

Theorem 4.1 |5l Theorem 3.1] Let u(x,y) be the potential density of a tran- 
sient Markov process with state space E. Then for every j3 > there exists 
a positive process {0 x ,x E E}, such that for every («i,...,a n ) in i?" and 
(xi, . . . , x n ) £ E, 



where U = {u(xi,Xj)}fj =1 . 

The proof of this theorem makes use of the fact that because u(x,y) is the 
potential density of a transient Markov process, U~ l is an M-matrix. It is a 
property of M-matrices that all its real eigenvalues are positive. Consequently 
all the real eigenvalues of U are positive. In addition the matrix rU r , see (|2.6|) . 
is a resolvent matrix and hence it has only non- negative entries. In this case it 
is trivial to see that U r is /3-positive definite for all /3 > 0. Therefore, Theorem 
14.11 follows from Proposition 12.11 

When a permanental process is determined by the potential of a transient 
Markov process we refer it it as an associated process. (This is ambiguous 
because the kernel of a permanental process is not unique.) So far the only 
permanental vectors and processes that we can show exist are associated pro- 
cesses and they have the additional property that they are infinitely divisible. 

There is another property of associated processes that enables us to sim- 
plify the statement of Theorem 13.11 In general we don't know whether the 
function d in this theorem is a metric on T. We show in the next lemma that 
when the permanental process is an associated process d is a metric on T. In 
this case we can add to the statement of Theorem 13.11 that the version 6' is 
continuous on (T, d) and that (|3. 19[) gives a uniform modulus of continuity on 



logF(S) 



(1/2) iog|/ + sr 



(4.2) 




(4.3) 



(T,d). 
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Lemma 4.2 Let {9 x ,x € E} be a (3-permanental process associated with a 
transient Markov process with potential u(x,y). Then 

d(x, y) = 4^/2/3 (u{x, x) + u(y, y) - 2 (u(x, y)u(y, x)) 1/2 ) 1/2 . (4.4) 
is a metric on E. 

The proof follows from the next two lemmas: 
Lemma 4.3 Let 

u ai fa \ 

a 2 v 71 . (4.5) 
Pi 72 w J 

where a\ot2 = a 2 , fiifo = b 2 and 7172 = c 2 . Suppose that A > 0. If B^ 1 is an 
M matrix then A~ l is an M matrix. 

Proof We first consider the case when abc > 0. We have 

^ / vw - c 2 /3 2 72 - oliw ai7i - (3 2 v 

B^ 1 = - — - /3i7i - a 2 w uw — b 2 a 2 f3 2 - 71 tt (4.6) 
deto loo 2 
«272 — Pl v «iPi — 72^ uv — a 





(u 


a 


M 




H 


a 


V 




1, .-( 




\ b 


c 


w J 





Since B 1 is an M matrix it's off diagonal terms are less than or equal to zero. 
This implies that 

gi := 7— > 1, 92 := — >!, 53 := -r > 1- (4.7) 
be ac ab 

We get this by considering the pairs B~^\, B^\ and B^l, B% \ and B^\, B%\ 
and using the hypothesis that B~ l is an M matrix. We now note that 

det*4 /uvw aw\ ^ A 6v\ ^ cu\ (4 8) 

a6c V a&c be J \ ac J V a&/ 

= (919293 - 9i) + (1 - 92) + (1 - 93) • 

One can check that on the range of (gi, 92, 93) allowed in (j4.7|) det A > unless 
two of the qi = 1, i = 1, 2, 3. (Set = 1 + ej, i = 1, 2, 3 and observe that (|4,8p 
is equal to e\e 2 + e 2 £3 + £1^3 + ^1^3-) 

Without loss of generality suppose q± = q 2 = 1. In this case, when abc > 0, 
B^ 1 can not be an M-matrix. We first note that if it is an M-matrix then 
{B~ 1 }i 2 = {-B _1 }2,i = 0. This holds since f5 2 r ) 2 < Oi\w and /?i7i < a 2 w. 
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Therefore, in order for b 2 c 2 = a 2 w 2 , which is the case when q\ = 1, we 
must have /3272 = ot\w and /3i7i = a 2 w. A similar argument shows that 
argument shows that {B~ 1 }i3 = {B^ 1 }^^ = 0. In addition {B~ 1 }i t i = 0, 
since q\ = q 2 = 1 implies that vw = c 2 . Consequently when q\ = q2 = 1 the 
first row of B^ 1 = 0. But this is impossible since an M-matrix is invertible. 

We have now established that under the hypotheses of this theorem, when 
abc > 0, det A > 0. 

Now suppose that abc = 0. It is easy to see from (|4.6p that this implies 
that either two of these three numbers must be zero or all three must be zero. 
For example, suppose that a\ = 0. Then by looking at B^\ we see that either 
72 or /?2 are equal to zero. Suppose ct\ = /32=0. Then 



det B = u(vw — c ) + ai/3i7i + 02/^272 — vfiifo — waia.2 (4.9) 

„2\ 



u[yw — c 2 

~2\ 



In addition, since in this case a = b = 0, det .4 = u(vw — c ). Therefore, 
since det B > we have det A > 0. A similar argument shows that when 
ot\ = 72 = and B~ 1 is an Af-matrix, we also have det A > 0. If all three 
numbers a, b, c are equal to zero, det A = uvw and so does det B, even though, 
once again, both B and B^ 1 can have non-zero off diagonal elements. In either 
case if B~ l is an M-matrix, det .A > 0. 
Since det A > 0, we have 

^ / vw — c 2 be — aw ac — bv \ 

A^ 1 = 7 be — aw uw — b 2 ab — cu I. (4.10) 

det A \ , , 2 / 

y ac-ow ao — cu uv — a ) 

If ^4 _1 is not an M matrix then one of the off diagonal terms of A^ 1 must be 
strictly positive. Without loss of generality we can suppose aw < be. Then 

aia 2 w 2 < /3i/3 2 7i72- (4.11) 

If this happens then one or both of the following inequalities must hold: 

ct\w < /3 2 72, a 2 w < /3i7i. (4.12) 

If either of these hold one of the off diagonal terms of B^ 1 must be positive. 
Therefore, B~ l is not an M matrix. This contradiction completes the proof. 
□ 
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Lemma 4.4 Let {6 x ,x £ R 1 } be a permanental process with kernel T(x,y). 
Let 

1 /2 

d(x, y) := (v(x, x) + T(y, y) - 2{T{x, y)T(y, x)) 1 / 2 ) . (4.13) 

Suppose that for any three different real numbers, x\, X2, x 3 , the 3x3 matrix 
r := {T(xi, Xj)}i t j=i t 2,3 has an inverse that is an M-matrix. Then d(x,y) is a 
metric on Br. 

Proof For any x±, X2, x 3 in R 1 it is easy to see that 

d(x 1 ,x 2 ) < d(x 1 ,x 3 ) + d(x 3 ,x 2 ). (4.14) 
This follows because by Lemma 14.31 the determinant of 

/ r( xi , xi ) (t(x 1 ,x 2 )t(x 2 ,x 1 )) 1 / 2 (r(xi,x 3 )r(x 3 ,xi)) 1 / 2 \ 

(T(x 1 ,X2)T(x 2 ,x 1 )) 1 / 2 T(X2,X 2 ) (T(x 3 ,x 2 )T(x2,x 3 )) 1 / 2 

V(r(x 1 ,x 3 )r(x 3 ,x 1 )) 1 /2 (T(x 3 ,X2)r(x 2 ,X 3 )) 1 / 2 T(X 3 ,X 3 ) J 

(4.15) 

is strictly positive. The determinant of the principle 2x2 minor is also strictly 
positive by \12\ (1.5)]. Therefore the matrix is positive definite and thus it 
is the covariance matrix of a Gaussian random variable which we denote by 
(G Xl ,G X2 ,G X3 ). Since 

d{x l ,x j ) = (E{G Xi -G X] ) 2 ) 1 ' 2 i,j = 1,2,3, (4.16) 
we get (|4.14p . □ 



Remark 4.2 Note that det£> — det.4 = (ai/3i7i — 02/3272) 2 so it is possible 
to have det£? > and det^4 < 0, but not when B~ l is an M-matrix. 

4.1 Continuity of associated permanental processes 

By Theorem 14.11 when u(x,y), x,y € T, is the potential density of a transient 
Markov process with state space T, a /3-permanental process 6/3 = {9 X ^ : x £ 
T}, with kernel n, exists for all f3 > 0. In addition, in this case, by Lemma [4. 21 
the function d in (|3.5[) is a metric on T. Using these facts we get the following 
extension of Theorem I3.lt 

Theorem 4.2 Let 9p = {0 x ^,x € T} be a /3-permanental process associated 
with a transient Markov process with potential u(x, y) satisfying sup^g^ u(x, x) < 
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oo. Let d be as given in j4-4\ ) an d ^ D denote the d diameter of T . Assume 
that there exists a probability measure on B(T, d) such that 

J S (D) < oo. (4.17) 

Then there exists a version 9' g = {9' x g , x € T} of Op such that for any xq € T 

|| sup ^||^ < Cp (we'Jl^ +supT(x,x)J 2 ~(D)) , (4.18) 
where Cp is a constant. If 

limJ^)=0, (4.19) 

there exists a version 0' g = {0' x g,x G T} of Bp that is continuous on (T,d). 
If ^4-19 ) holds and 



lim = oo, 4.20 

5^0 5 

then 

\9' - 9' | / \ 1 / 2 

lim sup -M y — < Cr ( sup 6' 8 ) a.s., (4.21) 

8^0 ^ yeT J^d(x, y)/2) W ,P J 

d(x,y)<8 

where Cr is a constant depending only on (3. 
Proof Given (3 choose p so that (3 < p/2. We have 

(p Stft - o yJl ) + (e^ p/2 _ B - e'^ p/2 _ B ) 4 £ (^ il/2 o) _ yil/2 o)) , (4.22) 

where 9p and Op are independent and where the notation (j) denotes inde- 
pendent copies. Let || • || be a pseudo norm on real valued functions on T. By 
(IQ2D . 

e\\ (e XiP - e VtP ) + (e' XtP/2 _ 8 - e' VjP/2 _p) \\ < p e\\9 x>1/2 - e y>1/2 \\. (4.23) 

Using the same argument we used in f)3.2T[) — (|3.3ip we get 

E\\9 x ,p - Ml < 2(p - P)E\\9 Xtl/2 - e yA/2 \\. (4.24) 
Similarly, we also note that starting with 

^ + ^ P /2^ = X>,i/2 ' ) (4.25) 
3=1 
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we can show that 



(4.26) 



Therefore, (I4TT81) follows from (f3T5l) . 
Using (|4.24|) we have 

E sup \0 XtP -9 Vi p\<2(p-P)E sup |^,i/ 2 -^,i/ 2 | (4.27) 

s,tET' _ s ' teT 

d(s,t)<6 d(s,t)<8 

for all finite subsets T' of T. By (|4.19p and (|3.17j) there exists a version 

e' p = {9' x ^xeT}oie^ 



YsmE sup \6' x p-6' y p\ 

s ,teT> 
d(s,t)<5 



(4.28) 



This implies that 8'^ is continuous on (T,d). 

By (|4.20p and (|3.19|) there exists a positive function g(5), satisfying lim^o 9($) 
0, such that 



' \f)' — f)' \ , \ 

\d(s,t)<5 



>31 sup^ 1/2 



< <?(£). (4.29) 



Using the argument in (I3.38h we get that there exists a constant Cp and a 
positive function g'(S), satisfying lim^o g'($) = 0, such that 



/ 



P 



\d(s,t)<5 



1 - 9' I / \l/2 



> C/3 sup9' 1/2 



<g'(S). 



(4.30) 



This implies that 

\0' — 9' I / X 1 / 2 

d(x,y)<5 

By |1 1|. Lemma 7.1.6] this is equivalent to (I4.2ip . 



(4.31) 



□ 
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5 Permanental processes associated with Levy pro- 
cesses 

It is interesting to see what kernels of permanental processes, that are not 
Gaussian squares, look like. We can do this by finding the kernels of perma- 
nental processes associated with Levy processes. Let X = {X tl t € R + } be a 
Levy process with characteristic function 

Ee iXX = e-^K (5.1) 

We refer to ip as the characteristic exponent of X. 

We assume that X has local times {L'f , (x,t) £ R x R+}. Set 

u To (x,y) = E* (L y To ) , (5.2) 

where To is the first hitting time of X at zero. We consider the transient 
Markov process X = {Xt} that is X killed at the first time it hits zero. The 
function UT (x,y) is the zero potential of X and thus is also the kernel of a 
permanental process. 



Lemma 5.1 



UT (x,y) = R(x,y) + H(x,y) and u To (y,x) = R(x,y)-H(x,y) (5.3) 
where 

1 f°° (1 - cos Ax - cos Ay + cosA(x - y))Keip{X) 

R{X > y) = *Jo W (5>4) 

and 

1 f°° (sin Ax - sin Ay - sinA(x - y))lmip(X) 

H(x,y) = --J o ^yp dX. (5.5) 

Proof It follows from [1 H Theorem 4.2.4], modified for non-symmetric Levy 
processes, that for x,y G R 1 

UT (x,y) = 4>{x) + 4>{-y) - 4>(x - y), (5.6) 
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where 



1 f 00 1 



1 f 00 (l — e iAx ) 

1 /- 00 (l - e iAa; ) V(-A) 1 f°° (l- e" iAx ) V(A) 



^ Jo l^(A)| 2 2^7 |V(A)| 

1 / p (i-e^)^(-A) ^ 



Since ?/>(A) = ip(—\) we see that 

* (x) = *U wop " A "/o W dA) - (5 - 8) 



Therefore, 



(x) + 1 /*°° (l-cosAx)^(A) ^ 

1 Z" 00 sin Xxlm ip(\) 



and 



2 Wo IVW 

We write 



■dA. (5.10) 



and 



Set 



/ \ _ U T (x, y) + u Tq (y, x) | u To (a, y) - u T() (y, x) 
u T {x,y) — - I - lo-J-J-J 



„\ _ u To (x,y) + u To (y,x) u To (x, y) - u To (y, x) 

u T Q {y,x) — - - . 



R(x,y) = UT ° ix > y)+ 2 UToiy > x) and H(x,y) = UT ^ y) . 

(5.13) 

Using (|5.6p and the following equations we get f)5.3j) (|5.5|) . □ 

We now consider the transient Markov process X = {Xt} that is X killed 
at £ a , an independent exponential time with mean a > . We assume that 
X has local times {Lf , (x, t) G R x R + } and set 

u a (x,y) = E*(Ly). (5.14) 
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The function u a (x,y) is the zero potential of X and thus is also the kernel of 
a permanental process. 

Lemma 5.2 For a > 0, the a-potential density 

u a (x,y) = R a (x,y) + H a (x,y) and u a (y, x) = R a (x, y) - H a (x, y), (5.15) 
where 

1 r°° co S \(x-y)Ke(a + TP(\)) 

Ra{x,y) = - : — 2 dX 5.16 



and 



= i /" "X'-^y^i fl, ( 5.i7) 



Proof By flSJJ 

/»oo 

e* A «p t (0,|/)dy = e-^W (5. 18) 



Therefore 



1 r 00 

(0,y) = — / e^e"*^ A\ (5.19) 
and, similarly to ()5.7[) 

u a (0 )2/ ) = / e - at p t (0,y)dt = — - dX (5.20) 

Jo 2vr y_ oc a + V(A) 



1 + ^)) 

vr \Jo |a + V(A)| 2 
1 f f 00 cosAy^e(a + V(A)) 



^ V7o |a + ^(A)| 2 

! ' 10 sin Ay Zm (a + ^(A)) 



/o |a + V(A)| 2 

Using the fact that u Q (x, y) = u a (0, y — x) we see that 

u a (x,y) + u a (y,x) 1 Z" 00 cos A(x — y)7£e (a + ^(A)) 



dX 



cos X(x-y)Ke (a + m)dX _ 
■A irj \a + ^{X)\ 2 

and 

u a (x,y) - u a (y,x) 1 f°° sin A(z - y) 1m (a + ip(X) 

o = _ / 1 — i 77TM2 dA ' 5 - 22 

which gives fj5 .151) . □ 
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Remark 5.1 Denote R and H in (j5.4j) and ()5.5|) by i?o an d -ffo- Since 



1— cos Ax — cos Ay+cos A(x — y) = (1— cos Ax)(l — cos Ay)+sin Ax sin Ay (5.23) 
and 

cos A(x — y) = cos Ax cos Xy + sin Ax sin Ay (5-24) 

we see that R a (x,y) is positive definite for all a > 0. Therefore they are the 
covariances of mean zero Gaussian processes G a = {G a (x),x € R 1 }. It is easy 
to see that for a > 0, G a is a stationary Gaussian process and that Go, has 
stationary increments, and that Go(0) = 0. Furthermore 

2 _2 f°° (1 - cosA(x - y))neip{X) 



E(G (x) - G (y)Y = - / ^ ^ dX (5.25) 

and when _E(G 2 (x)) = 1, 

E (G« W - G „ (s) )= = a r (1 ' cosA( r7,mf + ^ A)) v « > »• 

vr7 |a + V(A)| 2 

(5.26) 

Since R a (x,y) is positive definite for all a > 0, it follows from from (|5.3p that 
UT (x,y) and u a (x,y), a > 0, are positive definite. This is easy to see since 

UT (x,y) = Ro(x,y) +H (x,y) (5.27) 

and Hq(x,x) = 0, and Ho(x,y) = —Ho(y,x), and similarly for u a (x,y). 
Set 

d Ga (x,y) = (S(G a (x)-G Q (y)) 2 ) 1/2 a > 0. (5.28) 

We show in |12|. Lemma 5.5] that this metric is equivalent to d in (14. 4p . There- 
fore it follows from Theorem 14.21 an d |Hj Chapter 6], that the denominator 
in (|4.2ip for a permanental processes with kernel UT (x,y) or u a (x,y), a > 0, 
is the same, up to a constant multiple, as the uniform modulus of continuity 
of Gaussian processes with increments variance do a - 

5.1 FBMQ a/3 

There is a great deal of interest in fractional Brownian motion. Therefore, 
there should be considerable interest in permanental process for which the 
function R(x, y) in (I5.4p is the covariance of fractional Brownian motion. To 
study these processes we use Lemma 15. II to find the kernel of the permanental 
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process associated with the non-symmetric stable process of index 1 < a + 1 < 
2. This process has a Levy exponent of the form 

-0(A) = c|A| a+1 (l - i/3 sign (A) tan((a + l)vr/2)) (5.29) 

where c is an arbitrary constant which we set to 1, and |/3| < 1. We refer to 
this permanental process as FBMQ a,/3 . We use this notation because when 
P = 0, ^(A) is the Levy exponent of the Levy process associated with fractional 
Brownian motion of index a, (i.e. FBM). We add the Q, for quadratic, to 
denote the square of this process, as one does in the designation of the squared 
Bessel processes, (BESQ). 

Lemma 5.3 Let X = {Xt,t E R+} be a Levy process with characteristic given 
by 115. 29\) . Let UT (x,y) be as given in (5.ty) and write UT (x,y) and UT (y,x) 
as in 115. 3\) . In this case 

R(x, y) = C Q)/3 (\x\ a + \y\ a -\x- y\ a ) , (5.30) 

and 

H(x,y) 

= PC a fi ( sign (x)\x\ a - sign (y)\y\ a - sign (x - y)\x 

where 

-sin((q + l)f)r(-q) 
a3 tt(1 + /3 2 tan 2 ((a + 1)tt/2)) 

Consequently 

u To (x,y) = C a ^ ((1 - P sign (x)) \x\ a 

+ (1 - P sign (-y)) \y\ a - (1 - ft sign (x - yj) \x 

Proof Let c Qj(3 = (1 + f3 2 tan 2 ((a + l)vr/2))" 1 . Then 

(1 — cos\x)7Zeip(X) Cq,^ (1 — cos \x) 

WM 2 = |A| a+1 

and 

sinAxXmV'(A) c Qj/ g/3 sign (A) tan((a + l)ir/2) sin \x 

WM 2 = ; W+ l 



(5.31) 

-y\ a ), 

(5.32) 
(5.33) 

-y\ a )- 

(5.34) 
(5.35) 
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Therefore, by a change of variables 
sin XxTmip{X) 



o 



dX (5.36) 



2 



and 



s /„% f°° sin Ax „ 
= -c^ P tan((a + 1)tt/2) / d A 

= -c a ,/3/3tan((a + l)7r/2)|x| Q sign (x) y ~\^+T dX 
(1 -cosAx)^eV(A) / ,0O (l-cosA) 



o hKA)| 2 

It follows from [TJ 17.33.2] and an integration that 

r-oo i AX 



d\ = c a ^\x\ a (1 dX. (5.37) 



dA = - (sin ((a + + *cos ((a + T(-a). (5.38) 



Combining (|5.36p - (|5.38p we see that 

(1 — cos Ax) IZe ijj(X) 



-I 

n Jo 



mxw 

and 



d\ = C a ^ \x\ a (5.39) 



1 f°° smXxlmip(X) 

1 -dX = j3 sign (x)Cq, )/9 |x| . (5.40) 



Wo IV>(A)| 2 

Using (I5T39]1 and (pxlOD in ([53]) and (J53J) we get (OUjl and dOTTl . □ 

Remark 5.2 It follows from Remark 15.11 and [TT| (7.186)] that when ^ is 
FBMQ a,/3 then for any T < oo 



lim sup — — — W jp- —j- < Ct I sup 9 X ] a.s. (5.41) 

6^0 x , ye[lhT] (\x - y\ a logl/\x - y\) 1 ' 2 \xe[o,T] J 

d(x,y)<S 

In addition by |12[ Theorem 4.2] and Example 7.6.6, (1)] for xo ^ 

lim sup , J^fL%J < C fli/2 a . s . (5 .4 2) 

^dWoK^d^-^hloglogl/lx-xol) 1 ^ " 



1/2 



and 



lim sup t— — : — r ttt < C (5.43) 

^d^ox^ M a loglogl/|xj 
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A very interesting unanswered question is: can we find FBMQ" for 1 < 
a < 2 other than the square of a fractional Brownian motion itself? We give 
some examples to explain what we have in mind. Let G a = {G a (x),x € -R+} 
be fractional Brownian motion with index 1 < a < 2, i.e. G a (0) = and 

a 2 (x - y) := E(G a (x) - G a (y)) 2 = \x - y\ a . (5.44) 
It follows that the covariance 

f(x,y) := E(G a (x)G a (y)) = \ (\x\ a + \y\ a - (\x - y\ a )) . (5.45) 
Consider the kernel 

{f(x,y) + ea(x-y) x>y 
(5.46) 
T(x,y) - ea(x - y) x<y 

where e > is some small number. To simplify the notation set d(x, y) = 
d(x, y)/4y / 2/3. We can write 

d (x,y) = G(x,x) + G(y,y) -G{x,y) -g(y,x) (5.47) 
+G(x,y) + Q(y,x) - 2 (g(x, y)g(y, x)) 1/2 
= f(x, x) + f(y, y) - 2r(x, y) + {g l,2 {x, y) - g l ' 2 (y, x)) 2 

< o- 2 (x -y) + \g(x,y) -G(y,x)\ 

< (l + 2e)a 2 (x-y). 

Also we see from the second equality in (|5.47p that d (x,y) > a 2 {x — y). 

We do not know whether G(x, y) is the kernel of a permanental process. If 
it is, it would be a generalization of the square of fractional Brownian motion. 
It would be very useful to know whether such generalizations exist. 

Here is a very simple example of the above. Consider (15.451) when a = 1. 
This is the covariance of Brownian motion. Consider G(x,y) in this case for 
three values, < x < y < z. The kernel obtained is 




A e := \ x + elx-yl 1 ' 2 y y - e\x - y\ 1 ' 2 \. (5.48) 



When e = 0, A is the kernel of {B 2 ,B 2 ,B 2 } where {B x ,B y ,B z } are three 
values of standard Brownian motion. One can check that Aq 1 in an M-matrix. 
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This is not a surprise since we now that Brownian motion is an associated 
process. (See e.g. Lemma 2.5.1].) However, Aj 1 is not an M-matrix for 
any e ^ 0. It would be interesting to know whether A t can be the kernel of a 
permanental vector. 

Remark 5.3 One may consider kernels of the form (|5.46j) for any Gaussian 
process with stationary increments and try to answer the question, "Are they 
kernels of permanental processes?" 

6 Some problems 

1. Necessary and sufficient conditions for the continuity of perma- 
nental processes. This question could be phrased more specifically as, "Is 
Theorem 13.11 a necessary and sufficient condition for the continuity of 1/2- 
permanental processes?" It is clear that Theorem 13. II is best possible because 
it is a necessary and sufficient condition for the continuity of the square of a 
Gaussian process with continuous covariance; (see [I2j Section 7.2]). 

When one considers this question one is immediately aware of how many of 
the properties of Gaussian processes that are used in the proof of a necessary 
condition for continuity or boundedness are unknown for permanental pro- 
cesses. Let 9 = {9 x ,x € T} be a 1/2-permanental process. If sup xeT 9 x < oo, 
is E(sup xeT 9 X < oo)? What about || sup xgT OxW^ If the metric d in (13. 5h 
of one permanental process is uniformly less that the metric d of another per- 
manental process, does the continuity or boundedness of the process with the 
larger metric imply the continuity or boundedness of the process with the 
smaller metric? More generally, are there any analogues of Slepian's Lemma 
for permanental processes? Gaussian processes often satisfy many important 
zero one laws. What zero one laws do permanental processes satisfy? All of 
these are interesting questions. 

2. Independence of the components of a permanental vector. It 

follows from Lemma 13. II that if 9 := {9 x ,x € T} is a 1/2-permanental process 
with kernel T then for any pair x,y, 

cov(e x ,e y ) = 2T(x,y)r(y,x) (6.1) 

If T(x, y)T(y, x) = it follows from Lemma lBTTI ttiat 9 X and 9 y are independent. 
Another interesting question is, "If the elements of {9 x ,x S T} are pairwise 
independent are they independent?" 
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Suppose that T = {1, 2, 3} and the kernel of 9 is 




3.2) 



where abc ^ 0, and otherwise a, b, c are arbitrary. 

By the general formula for the moments of /3-permanental processes 

E \U 9 -A = Perm p{{^ x i)}l<iJ<n) (6-3) 

n 

= E acW n r (^a)) 

7T j = l 

where c(ir) is the number of cycles in the permutation tt of [l,n], (See 
[131 Proposition 4.2].) Therefore, if the kernel given by (16. 2\i defines a 1/2- 
permanental vector (8%, 82, 83), then we would have £(818283) =l/8+(l/2)abc. 
Whereas, if 81,82,83 were independent E(8\ 

8283) = 1/8. On the other hand by (16. ID 81,82,83 are pairwise independent. 

The kernel in (j6.2j) and its transpose are the only possibilities for 8\, 82,83 
to be pairwise independent but not independent. So in this case the problem 
becomes, can this be the kernel of a 1/2-permanental vector. We can show 
that it can not. But in we can not answer this question. In R^_ we can 
have kernels of the form 



/ 1 a \ 

1 b e 

1/ 

V d 1 J 



(6.4) 



3. Permanental vectors in R^_. Given a 3 x 3 matrix with components that 
satisfy (I3.2p - (l3.4p . is it the kernel of a permanental process? If the inverse of 
the matrix is an M-matrix, the answer is yes. At this point we don't know 
of any examples of matrices that are kernels of a permanental process with 
inverses that are not M-matrices, except, of course, when the permanental 
process is the square of a Gaussian process. We think that there are many 
examples. The class of associated vectors that are Gaussian squares is a very 
small subset of all vectors that are Gaussian squares. Roughly speaking, they 
don't include any vectors that are squares of Gaussian processes that are 
smoother than Brownian motion. If, again, roughly speaking, one can always 
find permanental vectors with kernels that deviate slightly from the kernels 



27 



of Gaussian squares that are infinitely divisible, one should also be able to 
find permanental vectors with kernels that deviate slightly from the kernels of 
Gaussian squares that are not infinitely divisible. (Here we are using infinitely 
divisible as a synonym for associated.) 
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